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The two-frequency problem of synchronization of the pulse train of a passively mode locked soliton 
laser to an externally injected pulse train is solved in the weak injection regime. The source and 
target frequency combs are distinguished by the spacing and offset frequency mismatches. Locking 
diagrams map the domain in the mismatch parameter space where stable locking of the combs is 
possible. We analyze the dependence of the locking behavior on the relative frequency and chirp 
of the source and target pulses, and the conditions where the relative offset frequency has to be 
actively stabilized. Locked steady states are characterized by a fixed source-target time and phase 
shifts that map the locking domain. 



I. INTRODUCTION 

Synchronization, or injection locking of a self-sustained 
oscillator to an external periodic signal is a paradigm of 
nonlinear dynamics [1 . The frequency of the target os- 
cillator adjusts itself to the external frequency for strong 
enough injection or small enough frequency mismatch, 
with a sharp threshold separating configurations with 
stable locked steady states from unlocked configuration. 
Injection locking in cw lasers has been extensively stud- 
ied experimentally and theoretically |2J . Injection locking 
is useful to stabilize the target laser versus phase diffu- 
sion and frequency drift, and can be achieved for very 
weak seeding provided the frequency mismatch is small 
enough. Strong injection leads to rich behavior, including 
bifurcations, multistability, excitability and chaos [5]. 

If the laser is mode locked, and the injected signal is 
pulsed, the target laser pulse train can injection- lock to 
the source pulse train. Unlike cw synchronization that 
implies entrainment of a single optical frequency, pulse 
train synchronization entails entrainment of two frequen- 
cies, namely the pulse repetition rate and the pulse phase 
shift per round trip. Since these two frequency determine 
the spacing and offset of the pulse train frequency comb, 
pulse train synchronization is equivalent to injection lock- 
ing of the target laser frequency comb to the comb of the 
source. In this way a very high-quality standard fre- 
quency comb source, that can be quite weak, could im- 
part its quality to the target laser, locking it to the stan- 
dard and reducing imperfections, such as pulse noise and 
frequency drift. A natural application is when the source 
is carrier-envelope phase locked, where injection locking 
locks the carrier-envelope phase of the target. When the 
repetition rates of the source and target are not matched 
but rationally related, synchronization results with a tar- 
get repetition rate that is the least common multiple of 
the source and free running target, facilitating generation 
of high-repetition rate frequency comb [3]. 

Compared to cw injection locking, there are few stud- 
ies of injection locking of mode locked lasers. The basic 



theory of the system was developed in [3] using soliton 
perturbation theory in the framework of the Haus mater 
equation model. Reference [3] derived estimates of the 
locking regime parameters, and studied numerically the 
pulse dynamics for some typical locked and non-locked 
parameters. The phenomenon was demonstrated experi- 
mentally in actively mode locked lasers in [5] and in pas- 
sively mode locked fiber lasers in (6j [7] . The effects of 
noise in injection locked mode locked lasers was studied 
in [8] . The combined effect of noise and injection in the 
mode locking transition was studied theoretically in [9 
and experimentally in [10_, showing critical phenomena 
with classical exponents. Rational injection locking was 
demonstrated experimentally in passively mode locked 
lasers in |3J [TT] and in harmonic mode locked lasers in 

nana. 

Here we focus on the injection locking threshold phe- 
nomenon. If the envelope of a pulse of the free-running 
target laser is ip[ Y (t), its pulse train waveform is 

^free-running = ^ 1pfr(t ~ nT r ) e my ° Tr (1) 
n 

where r r is the repetition rate, and vqt t is the pulse to 
pulse overall phase shift. Assuming that the source wave 
form is 

Source = J2 ^(t ~ «(l + u)r r ) e m ^)r r (2) 
n 

where tp so (t) is the source pulse wave form, the time slip 
rate u and phase slip rate v between the source and the 
target pulse trains measure the mismatch between the 
source and target frequency combs. Since the pulse to 
pulse phase shift is defined modulo 2ir, v is defined mod- 
ulo — , while the frequency u can take all real values. 

In standard single-frequency synchronization there is 
a threshold frequency mismatch below which the target 
is locked. At threshold the locked state loses stability 
in a saddle-node bifurcation. Here there are two mis- 
match parameters, u and v, so for a given ip so (t) there is 
a domain in the v, u plane of where the target frequency 
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comb is stably injection locked to the source comb. In 
the locking domain the target waveform is then 

V'injection-lockcd - £ lfa(t ~ "C 1 + U ) T r) j< V ° +V ^ (3) 
n 

where the ipn is the target pulse waveform in the presence 
of injection. 

The shape and area of the locking domain has to be 
calculated on a case-by-case basis. Here we study soliton 
mode locking where the free target pulse waveform is non- 
chirped soliton-like iph(t) — asech(crf), and the source 
pulse waveform is 

V-so(i) = Kasech{at)e l " t+ ^ (4) 

so that the injected pulse shape is an attenuated version 
of the free-running pulse shape; still, it follows from the 
analysis below that the locking properties depend weakly 
on the precise source pulse shape as long as it approxi- 
mately matches the target pulse. 

The parameters of the source pulse are the dimension- 
less injection strength k, frequency oj and chirp /?. The 
parameter lu is usually omitted from the definition of the 
pulse shape because it can be absorbed into the carrier 
frequency. Here it has to be spelled out explicitly because 
the difference between the peak gain frequencies of the 
source and target pulses is a crucial locking parameter. 

In most of our analysis we assume that the injection 
is weak, k« 1. As shown below, under this assumption 
the locked steady state and its stability depend only on 
the ratios - and - and the main effect of the injection 
on the target waveform is a shift A in its overall timing 
and a shift 5 in its overall phase relative to the source 
pulse, so that ip^ w asech(a(t — A))e lS . 

The main results of this work are the locking diagrams 
shown in figure [I] for several values of pairs io, j3. The 
boundaries of the locking domains are shown in black. In 
addition to the threshold, the locking diagrams also show 
the A and S level curves. Interestingly, although in gen- 
eral the locking domains are smaller for larger \ui\ and |/3| 
due to decreasing source-target pulse overlap, the shrink- 
ing is not uniform, and locking with some mismatch pa- 
rameters is possible only for nonzero relative frequency 
or chirp. When lj ^ there is a correlation between the 
timing and phase frequency mismatches, where for some 
v values there is a minimal u needed for synchronization 
and vice versa. 



II. THE INJECTION LOCKING MASTER 
EQUATION 

Our analysis is based on the Haus theory of anoma- 
lous dispersion passively mode-locked lasers with a fast 
saturable absorber [TJ], also called mode locked soliton 
lasers [15] . The Haus master equation is a multiple scale 
model of the electric field envelope ip(r,t), where t is the 
fast sub-roundtrip time scale, and r is the slow time scale 



describing the evolution of the wave form for an integer 
multiple of the roundtrip time r r . 

Using the gain band center as the reference frequency 
the equation of motion for the free target laser is 

d T i> = i(d? + 2H 2 )V +gbP] (i + id 2 ) + adVf)^ (5) 

where g is the overall gain, assumed to depend on the 
pulse energy J dt\ip\ 2 , 7 is the coefficient of parabolic 
spectral gain profile, and a is the transmissivity of the 
saturable absorber, including linear loss equal to — <r(0). 
The first two terms on the right-hand-side of equation 
(JE]) model the anomalous dispersion and Kerr nonlinear- 
ity, with coefficients fixed by appropriate choice of units. 
The conversion of the results to physical units and their 
implication for practical systems are discussed below. 

Since the Haus master equation describes the evolution 
of the pulse waveform between consecutive roundtrip pe- 
riods, the optical injection can be modeled by an inhomo- 
geneous seed term ipih(t, r), where ipih(t, nr r ) is the en- 
velope of the seed pulse waveform (in appropriate units) 
injected in the nth roundtrip in the frame defined by 
equation For the source waveform ^ the injection 
locking master equation is therefore 

d T yj = i(d? + 2\^\ 2 )^ + g(l +1 d^ 

+ a(\i^\ 2 )ij + Ao(t-UT)e^°+»)T 

(6) 

The injection locking master equation that was first de- 
rived in [?] is equivalent to equation The master 
equation was written in [4] in the frame of reference 
where the injected signal is exactly periodic, and there- 
fore contains an additional term in the master equation 
that arises from a frequency shift of the source waveform, 
that is accounted here by the parameter u>. 

The target waveform approaches an injection locked 
pulse train p]) if there is a stable solution of Eq. (JsJ) of 
the form tp\\(t — ur)e , ' t '" +l '' T . An unlocked solution typ- 
ically consists of two pulses, a strong pulse of free target 
waveform weakly perturbed by interaction with the seed, 
and a weak response waveform that is approximately a 
linearly filtered version of the seed pulse. The unlocked 
state is not stationary in any frame. 

III. SOLITON PERTURBATION THEORY 

In soliton mode locking we assume that the gain, loss 
and injection terms in the equation are small. Then the 
target laser waveform has a perturbed soliton component 

r) = a sech(a(t - S))e ip ^- S ^ (7) 

where 

S(t) = 2 J dr'p(r') + s (8) 

$(r)= [ T dT'(a(T') 2 -p(T') 2 )+ v (9) 



3 



In addition to the strong pulse, the solution of the master 
equation |6| includes a weak continuum component, as 
well as a response pulse if the target fails to lock to the 
source. We are interested in the soliton-like pulse com- 
ponent ([7]) whose properties govern the frequency comb 
locking. 

The slow time dynamics of the pulse parameters, am- 
plitude a, phase ip, frequency p, and timing s is obtained 
by projection of Eq. (|6| on the four NLS adjoint eigen- 
a, tp,p, s |16j that gives 



functions % , k 



d T a = 25(1 - I (a 2 + 3p 2 ))a + a (a) - % 



d T tp ■■ 

d T p ■ 

d T s - 



-\l9Va 2 + 14> S + A 

4Vv 



where 



ipk =Im / dtx k {t)^ so {t) 



a(a) = -Im / d^(i)a(|Vw(*)r)Vw(i) 



(10) 

(11) 
(12) 
(13) 

(14) 
(15) 



The seed and target are synchronized if the seed-target 
phase mismatch 5 = (z/ + v)t — $ and timing mismatch 
A = ut — S are locked. It follows that a locked state is 
a solution of Eqs. (10 13 1 with d T a 
and d T s 



d T p = 0, d T ip = 
The seed projections depend 



a" — p", and d T s = 2p. 
on the mismatches so that the injection problem reduces 
in this regime to the solution of four coupled nonlinear 
equations for the unknowns a,p, 5, A given the mismatch 
parameters u and v. 



The free-running version of equations ( 10 13 I describes 
mode locked soliton-like pulses with p — 0, a - 
a, is the solution of 



where 



G(a s ) = 2 fl (a,)(l 



|a 2 )a 



)a s +a{a s ) = (16) 



with free overall phase accumulation rate vq — a 2 s . The 
question of synchronization is most interesting when the 
injected signal is weak and the mismatches are small. 
Here weak injection means that the injected signal is 
weaker than the intrinsic target laser gain and loss 
effects (themselves weaker than the the dispersive ef- 
fects) so that the pulse width and frequency are only 
slightly shifted from their free running values, that is 
a = a — a s <C a s and p <C a s . In this case it is natural 
use units such that a s = 1, so that the locking equations 
become non-dimensional and linear in a and p 



d T a = —Aa — tp v 
d T S = v — 2a — ip a 
d T p = -Pp + $3 
d T A = u — 2p + ip p 



(17) 
(18) 
(19) 
(20) 



where A = —d a G(a s ), necessarily positive for pulse sta- 
bility, and P — ^jg(a s ). The smallness assumptions used 



IV. THE INJECTION-LOCKED STEADY STATE 

Our goal is to map the set of stable injection-locked 
steady states. In the steady state equations (17 20 1 can 



be further simplified since the direct forcing terms in the 
phase and timing equations of motion are smaller than 
the indirect forcing via the amplitude and frequency, as a 
consequence of a general property of soliton perturbation 
theory. Specifically, it follows from equations ( 17|19 1 
that |&| ~ I tpJJA ^> \ip a \, an d therefore the forcing term 
in equation ( 18 1 is negligible; the term ip p is negligible in 
equation (20iK)r analogous reasons. 



By the preceding arguments, given the injection wave- 
form and the timing and phase mismatches, a locked 
steady state is a solution (A, S) of the coupled nonlin- 
ear equations 



2k 



9k 

u=-f u (A,5) 



(21) 
(22) 



where 



Kf„ = Re(e iS J sech(t)i; so (t - A)dtj (23) 

nf u = Im(e a / sech(t) tanh(i)?/> so (i - A)di) (24) 

displaying explicitly the forcing strength coefficient k of 
the source waveform. The solution is stable if the eigen- 
values of the matrix (omitting zero elements for clarity) 



S = 



-P 



(25) 



V 



to derive Eqs. ( 17 20 1 imply that \-ip k \ < A, P < 1. 



have negative real parts. 

The results derived so far hold for weak injection with 
an arbitrary pulse shape. The precise locking proper- 
ties depend on the details of the pulse shape through the 
forcing functions f v , f u \ nevertheless, since the forcing 
functions are defined by overlap intervals, their values 
are of 0(1) for a source pulse shape that matches the 
free target pulse shape in width and frequency, and is 
not excessively chirped. We therefore choose the form 
Q, allowing for mismatch uj in center frequency and j3 
in chirp between source and target. Unlike the comb mis- 
match parameters v and u, u) and /3 are pulse parameters 
whose values are not fixed by the locking dynamics. We 
further comment below on the different roles played by 
these two pairs of parameters. 

The solution of the injection locking problem is pre- 
sented in figure [T] as a set of locking diagrams for different 
choices of ui and (3, that display the regions in parameter 
space where synchronization is stable. In all cases there is 
a stable locked state when v = u = 0, while for large v or 
u there are no stable locked states. Frequency comb lock- 
ing is made easier by stronger injection, and hindered by 
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FIG. 1: Frequency comb locking diagrams for an attenuated soliton source pulse envelope. The horizontal coordinate designates 
the mismatch v between the source and free target frequency comb offsets times the pulse amplitude restoring coefficient A 
divided by the forcing strength k . The vertical coordinate designates the mismatch u times the pulse frequency restoring 
coefficient P divided by k. All quantities are dimensionless, expressed in natural units as explained in the text. Points inside 
the locking domain, whose boundary is the black curve, designate values of mismatches that allow for stable injection locking, 
where the target comb is pulled and locked to the source comb. Curves of constant source-target time shift in the locked state 
are shown in green, and constant source-target phase shift curves are shown in orange. The level set values are indicated by 
hue, coded according to the color bars in the bottom-right part of the figure. The values of source-target frequency difference 
to and chirp /3 are (a) u> = 0., /3 = 0., (b) u = O.,0 = 0.5, (c) u = 1., = 0., (d) u = 1., /3 = 0.5,(e) ui = 3., /3 = 0., (f) w = 2., 
/3 = .2 



stronger restoring coefficients, so that the mismatches are 
multiplied by the restoring coefficients (v by A and u by 
P) and divided by k in the locking diagrams; the physical 
implications of the scaling are analyzed in the conclusions 



section below. The locked and non-locked regions in the 
v-u plane are separated by a curve where at least one the 
eigenvalues of S is zero, shown in black in the locking 
diagrams. It follows that at the locking boundary point 
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the Jacobian matrix ( g *y J is singular, so that at 

generic boundary points there is a saddle-node bifurca- 
tion, where stable locked solutions connect with unstable 
locked solutions on reaching the locking domain bound- 
ary from the interior. It also follows that the constant A 
and constant S curves, shown in Figs. [T|in green and or- 
ange hues (respectively) meet the boundary tangentially. 

In the non-chirped locking diagrams the A = 0, 5 = 
and A = 0, (5 = 7r are mapped to corners on the lock- 
ing domain boundary, since the Jacobian vanishes there 
identically. The corners are not structurally stable, and 
when the source is chirped one of the corners becomes 
smooth and the other a cusp. The boundary of a typical 
locking diagram is therefore smooth except for one cusp. 

The most evident conclusions from the locking dia- 
grams is that the locking region shrinks with both in- 
creasing frequency mismatch and increasing chirp. This 
results is to be expected since the overlap of the source 
and target pulses decreases when there is relative fre- 
quency or chirp, and the injection efficiency decreases. 
When uj and f3 become large the area of the locking do- 
main decreases fast to zero — as a power law for large P, 
and exponentially for large uj. Nevertheless, the shrink- 
ing is not monotonic, and some areas of the v-u plain are 
more susceptible to locking with moderate frequency mis- 
match and/or moderate chirp, see figure [2] This is made 
possible by adjustment of the locked phase shifts A and 
5 to compensate for greater mismatches. For non-zero uj 
the source waveform has a time-reversal odd component, 
so that the values of f u and f u become closer, and the 
admissible v and u intervals become correlated; the cor- 
relation becomes more pronounced for larger frequency 
mismatch. 

The locking diagrams display some discrete symme- 
tries. First, for any locked state v,w with shifts A, 5, 
there is an opposite solution —v, —w with shifts A, n + 5: 
at most one of these solutions is stable. For a time- 
reversal even source envelope the following identities hold 

f u (A, 6; uj, p) = /„(A, -S- -uj, -p) = 

f v (-A,S;-uj,p) = f v (-A,-S;uj,-p) (26) 
f w {A,5;uj,p) = -f w (A,-5;-u>,-P) = 

- f w (-A, 5; -uj, p) = f w (-A, -5; uj, (27) 

It follows that the locking diagrams for all ui,P values 
are determined by those with both non-negative: The 
locking diagram for ~ui,P is obtained from that of ui,P, 
by horizontal reflection and replacement of A, S by — A, 5, 
and the diagram for uj, — /3 is obtained from that of uj, P 
by a 7r rotation around the origin and replacement of A, S 
by — A, 7r — S. In particular, the locking boundaries of the 
P = diagrams are symmetric with respect to a half turn 
around the origin, and those of the uj — diagrams are 
symmetric with respect to horizontal reflection. 



K/P 




FIG. 2: Boundary of locking region for u> = 0., ft — 0. (black), 
uj — 0.5, P = 0. (green), and U) = 0., fi = 0.5 (blue). It is evi- 
dent that moderate source chirp and source-target frequency 
shift improves locking properties for some comb mismatch 
values. The axes are as in Fig. [l] 



V. CONCLUSIONS 

The main result of this work is a precise mapping of the 
locking of a passive mode locked laser frequency comb to 
that of an external pulse train. We focused on the weak 
injection regime where this effect is most striking. In con- 
trast with cw synchronization, frequency combs are spec- 
ified by two frequencies, that have to be locked together, 
making the locking problem two-dimensional. The dif- 
ference between the two combs was parametrized by the 
mismatch v between the comb offsets and the mismatch 
u between the comb spacings. The mismatches should be 
comparable or smaller to the forcing strength k divided 
by the restoring coefficients A and P (respectively) of the 
target laser dynamics. 

The locking diagrams show the synchronization behav- 
ior using dimensionless quantities: u is expressed as mul- 
tiples of a s , the inverse of the free target pulse width, and 
v, A, and P as multiples of the dispersion frequency — 
the dispersion coefficient P times a 2 s times the group 
velocity — that is typically much slower than a s . The 
injection-restoration ratios -j, p are a measure of the 
strength of the injection, here assumed small. For weak 
injection therefore, the repetition rate of the source and 
the target should be matched up to a fraction of the 
pulse bandwidth, and the offset frequencies have to be 
matched up to a fraction of the slow frequency unit. If 
k is increased, the locking domain becomes larger pro- 
portionately, until the forcing is no longer weak, when 
nonlinear effects not studied here may prevent further 
growth of the locking domain. 

Practically, the laser repetition rates can be tuned by 
changing the cavity length, and the pulse bandwidth 
makes the admissible u interval large enough for pas- 
sive control. On the other hand, the natural frequency 
scale of the offset mismatch v is much smaller, and the v 
tolerance is accordingly more limited. To understand the 
locking of pulse phase it should be recalled that v is actu- 
ally defined modulo where tr is the repetition rate, 
so that a physical v corresponds to a horizontal lattice in 
the locking diagrams. If v takes a random value between 
— -2- and ^r- (say) , there is a high probability of locking 
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if the lattice spacing ^ is of the order of or smaller than 
the range of offset mismatched in the diagram, and low 
probability otherwise. In fiber lasers with 0.5 ps pulses 
and standard fiber dispersion the slow frequency scale is 
about 1 MHz — much smaller than normal repetition rates 
of about 50 MHz. It is therefore difficult to obtain small 
enough v by passive control, and for this reason a sophis- 
ticated active stabilization of v was implemented in [3 
to achieve synchronization. 

The amplitude and frequency restoring coefficients A 
and P are hard to measure directly. A can be estimated 
by the relaxation rate of the gain medium, which can 
vary by many orders of magnitude between 10 2 sec -1 in 
Erbium and 10 9 sec -1 in diode lasers. P is equal to 
the product of the overall differential gain g and the 
gain bandwidth attenuation factor 7a 2 . g can be esti- 
mated from cw operation where it typically varies be- 
tween 10 6 sec _1 and 10 s sec -1 , while 7 is inversely pro- 
portional to the square of the bandwidth so that 7a 2 . 
varies from about 10 -2 for highly chirped pulses to 1 for 
bandwidth limited pulses. Weaker restoring coefficients 
allow for larger pulling of a and p from their free-running 
values, and therefore the v locking interval is larger for 
slow gain media and the u locking interval is larger for 
slower differential gain and larger gain bandwidth. 

An important conclusion from this discussion is that 
the locking properties depend strongly on the pulse du- 
ration and the group velocity dispersion. For fixed injec- 
tion efficiency k and restoring coefficients, the admissible 
v interval is proportional to /3 2 a^ and the admissible u 
interval is proportional to (3a s /j. It follows that injec- 
tion locking is particularly well-suited for high dispersion 
ultrafast lasers. 

Unlike the mismatch of comb parameters that acts as 
a barrier to synchronization, the matching of source and 
target pulse shapes affects the locking properties through 
the overlap integrals in equations ( 23 24 1 . It follows that 



the most relevant are the matching of the central fre- 
quency, chirp, and width of the pulses. Of these, the 



locking is most sensitive to frequency mismatch, since 
the overlap integrals drop sharply when the bands of the 
source and target do not match. A subtler effect is the 
deformation of the locking domain to an elongated shape, 
so that there is an 'easy' direction in the v-u plane where 
relatively large mismatches allow injection, which is not 
possible with comparable mismatches in other directions. 
In practice the frequencies are naturally matched if the 
two lasers use the same gain mechanism. 

The effect of relative chirp is to decrease the effective 
source bandwidth available for locking, so the area of the 
locking domain shrinks more slowly for large chirp than 
as for large frequency mismatch. The example of chirp 
also highlights the difference between the role of comb 
parameter and pulse parameter matching. Whereas for a 
given chirp there is always an area of the v-u plane that 
allows locking, no choice of pulse shapes lead to synchro- 
nization if the offset frequency mismatch v is significantly 
larger than k/A in natural units. These conclusions are 
in agreement with the experimental achievement of in- 
jection locking by tuning the chirp of the source pulse in 
®- 

We did not study in detail the effect of pulse-width 
mismatch between source target, but the effect is likely 
to be similar to that of chirp, as a gradual shrinking of 
the of the locking domain, but without the break- 
ing of the discrete symmetries. The study of the effect 
of noise which is liable to change some qualitative prop- 
erties [TJ Er 10J is beyond the scope of this work. The 
results derived here have natural implication in thresh- 
old phenomena in rational pulse injection locking and 
the resulting repetition rate multiplication that need to 
be studied separately. 

Acknowledgments J Reiner and E Olami participated 
in early stages of this project. The project was sup- 
ported by the Israel Science Foundation grant 1002/07, 
US-Israel Binational Science Foundation grant 2010442, 
D.K. was supported by an Australian Research Council 
Future Fellowship (FT110100513).. 



[1] Pikovsky, A, M Rosenblum, and J Kurths. "Phase Syn- 
chronization in Regular and Chaotic Systems: a Tuto- 
rial." Int J Bifurcat Chaos 10 2291-2306 (2000) 

[2] S. Wieczorek, B. Krauskopf, T. B. Simpson, and D. 
Lenstra, "The dynamical complexity of optically injected 
semiconductor lasers," Physics Reports, 416, 1-128, 
(2005). 

[3] D. Kielpinski and O. Gat, 'Phase-coherent repetition rate 
multiplication of a mode-locked laser from 40 MHz to 1 
GHz by injection locking,' Optics Express, 20, 2717-2724 
(2012) 

[4] M. Margalit, M. Orenstein, and H. Haus, "Injection lock- 
ing of a passively mode-locked laser," IEEE Journal of 
Quantum Electronics, 32, 155-160 (1996). 

[5] M. Margalit, M. Orenstein, G. Eisenstein, and V. 
Mikhaelshvili, "Injection locking of an actively mode- 
locked semiconductor laser," Optics letters, 19, 2125- 



2127, (1994). 

[6] M. Margalit and M. Orenstein, "Multimode effects on the 
evolution and long-term stability of a passively mode- 
locked laser under pulsed injection locking," IEEE Jour- 
nal of Quantum Electronics, 33, 710—718, (1998). 

[7] J. Dorring, L. Yan, D. A. Satorius, and Y. J. Chen, "Injec- 
tion locking of an additive-pulse mode-locked fiber laser," 
Photonics Technology Letters, 14, 1497—1499 (2002). 

[8] M. Margalit, M. Orenstein, and H. Haus, "Noise in pulsed 
injection locking of a passively modelocked laser," IEEE 
J Quantum Elect, 32, 796—801 (1996). 

[9] Weill R, A Rosen, A Gordon, O Gat, B Fischer, 'Critical 
Behavior of Light', Phys. Rev. Lett., 95, 013903 (2005). 
[10] A. Rosen, R. Weill, B. Levit, V. Smulakovsky, A. Bekker, 
and B. Fischer, "Experimental Observation of Critical 
Phenomena in a Laser Light System," Phys. Rev. Lett., 
105, 013905, (2010). 



7 



[11] M. Margalit, M. Orenstein, and G. Eisenstein, "High- 
repetition-rate mode-locked Er-doped fiber lasers by har- 
monic injection locking," Optics letters, 20, 1791-1793, 
(1995). 

[12] X. Wang, H. Yokoyama, and T. Shimizu, "Synchro- 
nized harmonic frequency mode-locking with laser diodes 
through optical pulse train injection," IEEE Photonics 
Technology Letters, 8, 617-619, (1996). 

[13] H. Kurita, T. Shimizu, and H. Yokoyama, IEEE Journal 



of Selected Topics in Quantum Electronics, 2, 508-513, 
(1996). 

[14] H. Haus, "Mode-locking of lasers," Ieee J Sel Top Quant, 
6, 1173-1185 (2000). 

[15] JN Kutz, "Mode-Locked Soliton Lasers," SIAM review, 
48, 629, (2006). 

[16] H. A. Haus and W. S. Wong, "Solitons in optical com- 
munications," Rev Mod Phys, 68, 423-444 (1996). 



